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Abstract 

The evolution of spherically symmetric unstable scalar field configura- 
tions ("bubbles") is examined for both symmetric (SDWP) and asymmet- 
ric (ADWP) double-well potentials. Bubbles with initial static energies 
Eq < Ecrit, where i^^rit is some critical value, shrink in a time scale deter- 
mined by their linear dimension, or "radius" . Bubbles with Eq > E^nt evolve 
into time-dependent, localized configurations which are very long-lived com- 
pared to characteristic time-scales in the models examined. The stability 
of these configurations is investigated and possible applications are briefiy 
discussed. 
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A remarkable consequence of nonlinear field theories is the existence of localized, 
non-singular solutions of the classical equations of motion which are non-dissipative. In 
general, these solutions can be time-dependent or static. As is well-known, the existence 
and simplicity of static solutions is severely constrained by dimensionality [jl|. For a 
single self-interacting real scalar field 0(x, t) such solutions are only possible in (1+1)- 
dimensions. More realistic (3+l)-dimensional static solutions must invoke more than one 
field, as in the case of the 't Hooft-Polyakov monopole [||. 

Given their relevance to the study of nonperturbative effects in field theories, static, 
non-dissipative solutions have been, for the last twenty years or so, the focus of most 
efforts in the study of nonlinear solutions in classical field theories. However, within the 
last decade, the possibility that spontaneous symmetry breaking occurred in the early 
Universe has called for a better understanding of time-dependent phenomena in the 
context of relativistic field theories. For example, the dynamics of cosmological phase 
transitions naturally invokes out-of-equilibrium conditions, with fields interacting with 
themselves and with a hot plasma in the background of an expanding Universe @]. 

In the present work, the possibility that time- dependent, localized, non-dissipative so- 
lutions exist in the context of simple (3+l)-dimensional scalar field theories is examined. 
In particular, the focus will be on models involving only a single real scalar field with 
self-interactions dictated by a double-well potential. Since it is known that for a sym- 
metric double-well potential (SDWP) all field configurations are unstable, it is possible 
to obtain the lifetime of a given spherically symmetric field configuration by numerically 
evolving the equation of motion. By addopting this procedure, it has been shown in 
the mid-seventies that certain configurations evolved into a state which was considerably 
long-lived (even though the lifetimes measured then were not very accurate). These time- 
dependent solutions were called "pulsons" However, not much has been done in order 
to further explore the properties of pulsons. A few exceptions, which are mainly related 
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to the existence of these solutions for the Sine-Gordon potential, different symmetries, 
and somewhat contrived stability studies, are listed in Ref. ^. In fact, recent studies of 
bubble evolution in SDWPs overlooked the existence of pulsons [0. Also, their existence 
has never been investigated for ADWPs. 

Here, it will be argued that the existence of pulsons is a very general feature of models 
with both symmetric and asymmetric potentials, depending only on the initial amplitude 
and energy of the configuration. It will always be assumed that the initial configuration 
interpolates between the two minima of the potential (bubbles with amplitude below a 
certain value cannot evolve into pulsons, as will be clear later), and that it is smooth 
enough, say, 0(r, t = 0) ~ exp[— r^/i^g], or ~ tanh(r — _Ro)) with Rq the initial "radius". 
(From now on these configurations will be called the gaussian and the tanh-bubbles.) 
It is remarkable that during the nonlinear evolution of these configurations a regime of 
dynamical stability may be achieved in which energy is practically conserved within a 
localized region, due to some as yet unclear conspiracy between the many degrees of 
freedom involved. 

The action for a real scalar field is, S'[<^] = / d'^x \d^(j)d^(j) — Vs{K){,(t)) , where the 
subscripts S and A stand for the SDWP and ADWPs, with, 1/s(0) = i (^^ - ^ )^ and 
Ka(0) = — ^mcfr' + ^(f)'^. A field configuration 0o(x, t) is a solution of the equation 
of motion (an overdot denotes partial time derivative), — = — and has an 
energy, 



E Unl = d'x 



For a sufficiently well-localized spherically symmetric configuration (j)Q{r,t) (the case 
of interest here) with linear "size" ~ Rq, the integral over all space can be restricted to 
a spherical volume containing the configuration, 

E [0o] (t) = EK{t) + Es{t) + Ev{t) , (2) 
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where the kinetic, surface, and volume energies are defined respectively by (e > 2) 

reRo . reRo „ r^Ro 

EKit) = 27T dr 0^, Esit) = 2n dr , Ey(t) = An dr V{<Po) • 

JO Jo Jo 

(3) 

Introducing the dimensionless variables $(r, t) = ^/Xcp^r, t) /m, p = m r, and r = mt, 
the equation of motion becomes, for both SDWP (and ADWP), 

where a = a/\/\. Note that for the ADWP, a solution where $ = is the local 
minimum is possible only if > 9/2. The equation above was solved numerically using 
a finite differencing method fourth order accurate in space and second order accurate in 
time. Since the problem is two-dimensional a very fine grid could be used. By taking 
the lattice spacing to he h = 10~^ and the time step to be = 5 x 10~^, energy was 
conserved throughout the evolution to one part in 10^. 

Let us concentrate on the SDWP for now. In this case all bubbles are unstable, since 
there is no gain in volume energy in going from one vacuum to another. Take the vacuum 
to be at $ = —1 and consider configurations which interpolate between the two minima. 
Both thick and thin-wall bubbles will be considered. In Fig. 1 the energy within a 
spherical shell, Eq. ||, is shown for initial configurations, ^o{p,0) = — tanh(p — po), with 
po = 3 and po = 10, and for ^o{p,0) = 2exp [— p^/Po] — 1, with po = 4 and po = 8. 
Note the existence of an extended period of stability, of duration ~ 10^ — lO^m^^, where 
practically no energy is radiated away. Thus, a necessary condition (which was verified 
numerically) for the existence of this pseudo-stable behavior is, from Eq. |^ 

dEK^_(dEs dEy\ 
dt \ dt dt ) ' ^ ' 

An interesting point is that the value of the energy at the plateau is fairly independent 
of the initial configuration. This suggests the existence of an attractor-like configuration 



in field space which is approached in the course of the bubble's evolution. An extensive 
(but not exaustive) search indicates that only for configurations of initial energies Eq > 
60m/ A = Ecrit this behavior occurs. This corresponds to a gaussian bubble of radius 
Rq ~ 2.4/m. For smaller energies, the evolution is well-fitted by the relation E{t) = 
Eoexp[—t/TL{Ro)], where tl{Ro) is the bubble's lifetime. For gaussian bubbles of radii 
rriRQ = 1, 2, mTL^Ro) = 3.5, 12, respectively. 

From Fig. 1 it is clear that the evolution of large enough bubbles can be divided into 
three stages. First the bubble sheds its initial energy by quickly shrinking into a thick- 
wall bubble of energy roughly E ~ 50m/ X. For thin wall bubbles, like the tanh-bubble 
with po ^ 1) this shrinking is well described by the relativistic motion of the bubble 
wall. (See, e.g., Ref. |^.) The field then settles into the pulson configuration. At any 
given time, the configuration is well-approximated by a half-gaussian, but with softer 
asymptotic behavior, $(p oo,r) ~ exp(— p). The field is localized within a small 
volume with linear size ~ 3/m, while its amplitude is rapidly oscillating about a value 
between the two minima. (In this author's opinion, a better name for these configurations 
would be "oscillons".) Finally, during the last stage of the evolution, the amplitude of 
oscillations decreases and the pulson quickly radiates its remaining energy away. 

In order to gain some insight into the properties of the pulson consider the behavior 
of radial perturbations about $o(p, r), defined as $(p, r) = $o(p, r) + (5$(p, r). Since 
^o{p, t) satisfies the equation of motion, expanding (5$(p, r) in normal modes, (5$(p, r) = 
Re J2n 0n(p)exp(zti;„r), it is found that the amplitudes 0„(p) satisfy the radial Schrodinger 
equation, (only the i = mode will be considered here), 

1— + V{p, ro)0„ = Xn(pn , (6) 

dp"^ p dp 

where V{p,tq) = 3$q(p, tq) — 3, and A„ = cj^ — 2. A„ is introduced to make sure that 
V{p oo) — > 0. The stability of a given configuration at time Tq, ^o{p, Tq), is determined 
by the lowest eigenvalue being positive, cUq > 0, or Aq > — 2. 
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Clearly the general problem is quite complicated, as the potential V{p, r) is time- 
dependent. A reasonable simplification is to solve this equation for a given time tq for 
which the field configuration ^q{p,tq) is obtained by evolving the equation of motion. 
This implicitly assumes that ^o{p,t) varies slower than 5$(p, r). Otherwise, unstable 
modes would not have enough time to grow and destabilize the pulson. One can then find 
the eigenvalue for a succession of snapshots and thus investigate its time evolution. The 
Schrodinger equation was solved using the shooting method P]. To make sure the method 
worked, the equation was solved for the Coulomb potential and for the "kink potential" 
V{p) = 3cosh-^(p/V2), since in both cases the eigenvalues are known analytically 0. 
(For the kink case, one must recall that the first eigenvalue in 3d corresponds to the 
second in Id, Eo{3d) = Ei{ld) = —1/2, due to the boundary conditions at the origin.) 

The results are shown in Fig. 2 during part of the pulson stage. Clearly, the pulson 
is stable against small radial perturbations. (No instability was detected during the 
whole pulson stage using the above method. An argument to explain the pulson's final 
disappearance is advanced shortly.) 

Next, the existence of pseudo-stable behavior for ADWPs is investigated. The asym- 
metry is controlled by the dimensionless coupling a. As two examples, consider the 
nearly degenerate case, a = 2.16, and the non-degerate case, a = 2.23. In the context 
of first-order phase transitions (finite temperature), a state initially localized at $ = 
will decay by the nucleation of bubbles larger than the critical bubble, which is an ex- 
tremum of the 0(3)-invariant Euclidean action. The critical bubble is thus the energy 
barrier for vacuum decay. For a given value of a it is easy to obtain the critical bubble 
and its energy numerically, by using the shooting method. (Now, the method finds the 
value of the field in the bubble's core.) For the two values of the asymmetry above, the 
energies and radii of the critical bubble are, respectively, E{a = 2.16) ~ 6.12 x lO^m/A, 
R{a = 2.16) ~ 19.6/m, and E{a = 2.23) ~ 1.17 x lO^m/A, R{a = 2.23) ~ 7.9/m. The 
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interest here is in investigating the evolution of sub-critical bubbles, to see if they display 
the pseudo-stable behavior observed for the SDWP. The results are shown in Fig. 3 
where the energy within a spherical shell surrounding the initial configurations is shown 
as a function of time. For each value of a, two tanh-bubbles with initial radii po = 3 and 
po = 4 were examined. Again, the existence of very long-lived pulsons is observed, with 
Ecrit and lifetimes depending on the asymmetry. In Fig. 4 the phase-space portrait of the 
pulson's core (p = 0), for a = 2.16 and po = 4 is displayed for r > 1000. Note the simi- 
larity with the motion of a "damped" (an)harmonic oscillator. During the pulson stage, 
the motion is restricted to a band in phase-space. As the pulson becomes unstable, it 
will spiral around the final stable point $ = and $ = 0; As energy is gradually radiated 
away, the maximum possible amplitude (that is for $ ~ 0) is driven below a critical value 
for stability. Bubbles with initial amplitudes below this value, roughly about the maxi- 
mum of V{^), and small kinetic energy will quickly shrink. Note that this is the value 
below which the potential is approximately parabolic, so that nonlinear effects become 
subdominant. (Of course, small amplitude but high velocity bubbles can still escape the 
attractive well centered at $ = 0, as during the pseudo-stable pulson stage.) 

The present work raises many questions for future investigation. Apart from a more 
detailed study of the pulson's properties and stability, it is still not clear why the initial 
configurations evolve into the pulson stage. The pulson's behavior is suggestive of some 
sort of nonlinear resonance effect occurring between the different modes of the field. It 
should be possible to separate the field into short and long wavelengths, with a cutoff 
around m~^. The short wavelength modes would act as a perturbation, which due to 
the nonlinear coupling might induce the observed behavior. Apart from their potential 
interest for, among other topics, nonlinear optics and long Josephson junctions [p!0| , 
these solutions are also of interest in the context of cosmological phase transitions [Q. 
Consider two possible regimes of interest, defined by the ratio of time scales TZ = t^/tjj ^ 
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{(j)/(J)) /{rripi/T'^) ~ T/rripi, where T is the temperature and mpi is the Planck mass. For 
TZ ~ 10^^, corresponding to the GUT scale, the extended lifetime of one sub-critical 
bubble can be relevant to the dynamics of a first-order transition. The Universe cools off 
quite fast and if the sub-critical bubble lives long enough it can become a critical bubble. 
For TZ ~ 10^^^, corresponding to the electroweak scale, many sub-critical bubbles can 
be present before a critical bubble nucleates |TT[. Given their long lifetimes, they could 
serve as nucleation sites for critical bubbles, very much like impurities in ordinary phase 
transitions, speeding the process of vacuum decay considerably. Finally, in the context 
of late-time phase transitions, it is desirable to have bubbles live longer in order for 



matter accretion to be efficient [|12|. These and related questions are currently under 
investigation. 
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Figure 1. The energy within a spherical shell surronding several initial configurations as 
a function of time for the SDWP. Larger bubbles have shorter lifetimes. 

Figure 2. Time evolution of lowest eigenvalue. For clarity, points with uJq > 2 were dis- 
played as a;Q = 2. 

Figure 3. The energy within a spherical shell surronding several initial configurations as 
a function of time for the ADWP. 

Figure 4. Phase-space evolution of pulson's core for a — 2.16, po — 4, and r > 1000. 
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